Evolution of genes involved in the unusual genitals of the bear macaque, Macaca arctoides.

Genital divergence is thought to contribute to reproductive barriers by establishing a "lock-and-key" mechanism for reproductive compatibility. One such example, Macaca arctoides, the bear macaque, has compensatory changes in both male and female genital morphology as compared to close relatives. M. arctoides also has a complex evolutionary history, having extensive introgression between the fascicularis and sinica macaque species groups. Here, phylogenetic relationships were analyzed via whole-genome sequences from five species, including M. arctoides, and two species each from the putative parental species groups. This analysis revealed ~3x more genomic regions supported placement in the sinica species group as compared to the fascicularis species group. Additionally, introgression analysis of the M. arctoides genome revealed it is a mosaic of recent polymorphisms shared with both species groups. To examine the evolution of their unique genital morphology further, the prevalence of candidate genes involved in genital morphology was compared against genome-wide outliers in various population genetic metrics of diversity, divergence, introgression, and selection, while accounting for background variation in recombination rate. This analysis identified 67 outlier genes, including several genes that influence baculum morphology in mice, which were of interest since the bear macaque has the longest primate baculum. The mean of four of the seven population genetic metrics was statistically different in the candidate genes as compared to the rest of the genome, suggesting that genes involved in genital morphology have increased divergence and decreased diversity beyond expectations. These results highlight specific genes that may have played a role in shaping the unique genital morphology in the bear macaque.

Altered fire regimes modify lizard communities in globally endangered Araucaria forests of the southern Andes.

Wildfire regimes are being altered in ecosystems worldwide. The density of reptiles responds to fires and changes to habitat structure. Some of the most vulnerable ecosystems to human-increased fire frequency are old-growth Araucaria araucana forests of the southern Andes. We investigated the effects of wildfires on the density and richness of a lizard community in these ecosystems, considering fire frequency and elapsed time since last fire. During the 2018/2019 southern summer season, we conducted 71 distance sampling transects to detect lizards in Araucaria forests of Chile in four fire "treatments": (1) unburned control, (2) long-term recovery, (3) short-term recovery, and (4) burned twice. We detected 713 lizards from 7 species. We found that the density and richness of lizards are impacted by wildfire frequency and time of recovery, mediated by the modification of habitat structure. The lizard community varied from a dominant arboreal species (L. pictus) in unburned and long-recovered stands, to a combination of ground-dwelling species (L. lemniscatus and L. araucaniensis) in areas affected by two fires. Araucaria forests provided key habitat features to forest reptiles after fires, but the persistence of these old-growth forests and associated biodiversity may be threatened given the increase in fire frequency.

Mitonuclear interactions and introgression genomics of macaque monkeys (Macaca) highlight the influence of behaviour on genome evolution.

In most macaques, females are philopatric and males migrate from their natal ranges, which results in pronounced divergence of mitochondrial genomes within and among species. We therefore predicted that some nuclear genes would have to acquire compensatory mutations to preserve compatibility with diverged interaction partners from the mitochondria. We additionally expected that these sex-differences would have distinctive effects on gene flow in the X and autosomes. Using new genomic data from 29 individuals from eight species of Southeast Asian macaque, we identified evidence of natural selection associated with mitonuclear interactions, including extreme outliers of interspecies differentiation and metrics of positive selection, low intraspecies polymorphism and atypically long runs of homozygosity associated with nuclear-encoded genes that interact with mitochondria-encoded genes. In one individual with introgressed mitochondria, we detected a small but significant enrichment of autosomal introgression blocks from the source species of her mitochondria that contained genes which interact with mitochondria-encoded loci. Our analyses also demonstrate that sex-specific demography sculpts genetic exchange across multiple species boundaries. These findings show that behaviour can have profound but indirect effects on genome evolution by influencing how interacting components of different genomic compartments (mitochondria, the autosomes and the sex chromosomes) move through time and space.

Mitogenomics of macaques (Macaca) across Wallace's Line in the context of modern human dispersals.

Wallace's Line demarcates a biogeographical boundary between the Indomalaya and Australasian ecoregions. Most placental mammalian genera, for example, occur to the west of this line, whereas most marsupial genera occur to the east. However, macaque monkeys are unusual because they naturally occur on both western and eastern sides. To further explore this anomalous distribution, we analyzed 222 mitochondrial genomes from ∼20 macaque species, including new genomes from 60 specimens. These comprise a population sampling of most Sulawesi macaques, Macaca fascicularis (long-tailed macaques) specimens that were collected by Alfred R. Wallace and specimens that were recovered during archaeological excavations at Liang Bua, a cave on the Indonesian island of Flores. In M. fascicularis, three mitochondrial lineages span the southernmost portion of Wallace's Line between Bali and Lombok, and divergences within these lineages are contemporaneous with, and possibly mediated by, past dispersals of modern human populations. Near the central portion of Wallace's Line between Borneo and Sulawesi, a more ancient dispersal of macaques from mainland Asia to Sulawesi preceded modern human colonization, which was followed by rapid dispersal of matrilines and was subsequently influenced by recent interspecies hybridization. In contrast to previous studies, we find no strong signal of recombination in most macaque mitochondrial genomes. These findings further characterize macaque evolution before and after modern human dispersal throughout Southeast Asia and point to possible effects on biodiversity of ancient human cultural diasporas.

Correction to 'Speciation over the edge: gene flow among non-human primate species across a formidable biogeographic barrier'.

[This corrects the article DOI: 10.1098/rsos.170351.].

Speciation over the edge: gene flow among non-human primate species across a formidable biogeographic barrier.

Many genera of terrestrial vertebrates diversified exclusively on one or the other side of Wallace's Line, which lies between Borneo and Sulawesi islands in Southeast Asia, and demarcates one of the sharpest biogeographic transition zones in the world. Macaque monkeys are unusual among vertebrate genera in that they are distributed on both sides of Wallace's Line, raising the question of whether dispersal across this barrier was an evolutionary one-off or a more protracted exchange-and if the latter, what were the genomic consequences. To explore the nature of speciation over the edge of this biogeographic divide, we used genomic data to test for evidence of gene flow between macaque species across Wallace's Line after macaques colonized Sulawesi. We recovered evidence of post-colonization gene flow, most prominently on the X chromosome. These results are consistent with the proposal that gene flow is a pervasive component of speciation-even when barriers to gene flow seem almost insurmountable.

Genetic variants related to disease susceptibility and immunotolerance in the Duffy antigen receptor for chemokines (DARC, Fy) gene in the black lion tamarin (Leontopithecus chrysopygus, primates).

The DARC (Duffy antigen receptor for chemokines) gene encodes the DARC protein, which serves multiple roles in the immune system, as a binding site for the malarial parasites Plasmodium vivax and Plasmodium knowlesi, a promiscuous chemokine receptor and a blood group antigen. Variation in DARC may play particularly significant roles in innate immunity, immunotolerance and pathogen entry in callitrichines, such as the black lion tamarin (Leontopithecus chrysopygus). We compared amino acid sequences of DARC in the black lion tamarin (BLT) to non-human Haplorhine primates and Homo sapiens. Consistent with prior studies in other Haplorhines, we observed that the chemokine receptor experiences two opposing selection forces: (1) positive selection on the Plasmodium binding site and (2) purifying selection. We observed also that D21N, F22L, and V25L differentiated BLT from humans at a critical site for P. vivax and P. knowlesi binding. One amino acid residue, F22L, was subject to both positive selection and fixation in New World monkeys, suggesting a beneficial role as an adaptive barrier to Plasmodium entry. Unlike in humans, we observed no variation in DARC among BLTs, suggesting that the protein does not play a role in immunotolerance. In addition, lion tamarins differed from humans at the blood compatibility Fya /Fyb antigen-binding site 44, as well as at the putative destabilizing residues A61, T68, A187, and L215, further supporting a difference in the functional role of DARC in these primates compared with humans. Further research is needed to determine whether changes in the Plasmodium and Fya /Fyb antigen-binding sites disrupt DARC function in callitrichines.

Multicopy gene family evolution on primate Y chromosomes.

BACKGROUND: The primate Y chromosome is distinguished by a lack of inter-chromosomal recombination along most of its length, extensive gene loss, and a prevalence of repetitive elements. A group of genes on the male-specific portion of the Y chromosome known as the "ampliconic genes" are present in multiple copies that are sometimes part of palindromes, and that undergo a form of intra-chromosomal recombination called gene conversion, wherein the nucleotides of one copy are homogenized by those of another. With the aim of further understanding gene family evolution of these genes, we collected nucleotide sequence and gene copy number information for several species of papionin monkey. We then tested for evidence of gene conversion, and developed a novel statistical framework to evaluate alternative models of gene family evolution using our data combined with other information from a human, a chimpanzee, and a rhesus macaque. RESULTS: Our results (i) recovered evidence for several novel examples of gene conversion in papionin monkeys and indicate that (ii) ampliconic gene families evolve faster than autosomal gene families and than single-copy genes on the Y chromosome and that (iii) Y-linked singleton and autosomal gene families evolved faster in humans and chimps than they do in the other Old World Monkey lineages we studied. CONCLUSIONS: Rapid evolution of ampliconic genes cannot be attributed solely to residence on the Y chromosome, nor to variation between primate lineages in the rate of gene family evolution. Instead other factors, such as natural selection and gene conversion, appear to play a role in driving temporal and genomic evolutionary heterogeneity in primate gene families.

Patterns of genetic variation and the role of selection in HTR1A and HTR1B in macaques (Macaca).

BACKGROUND: Research has increasingly highlighted the role of serotonin in behavior. However, few researchers have examined serotonin in an evolutionary context, although such research could provide insight into the evolution of important behaviors. The genus Macaca represents a useful model to address this, as this genus shows a wide range of behavioral variation. In addition, many genetic features of the macaque serotonin system are similar to those of humans, and as common models in biomedical research, knowledge of the genetic variation and evolution of serotonin functioning in macaques are particularly relevant for studies of human evolution. Here, we examine the role of selection in the macaque serotonin system by comparing patterns of genetic variation for two genes that code for two types of serotonin receptors - HTR1A and HTR1B - across five species of macaques. RESULTS: The pattern of variation is significantly different for HTR1A compared to HTR1B. Specifically, there is an increase in between-species variation compared to within-species variation for HTR1A. Phylogenetic analyses indicate that portions of HTR1A show an elevated level of nonsynonymous substitutions. Together these analyses are indicative of positive selection acting on HTR1A, but not HTR1B. Furthermore, the haplotype network for HTR1A is inconsistent with the species tree, potentially due to both deep coalescence and selection. CONCLUSIONS: The results of this study indicate distinct evolutionary histories for HTR1A and HTR1B, with HTR1A showing evidence of selection and a high level of divergence among species, a factor which may have an impact on biomedical research that uses these species as models. The wide genetic variation of HTR1A may also explain some of the species differences in behavior, although further studies on the phenotypic effect of the sequenced polymorphisms are needed to confirm this.

Reduced representation genome sequencing suggests low diversity on the sex chromosomes of tonkean macaque monkeys.

In species with separate sexes, social systems can differ in the relative variances of male versus female reproductive success. Papionin monkeys (macaques, mangabeys, mandrills, drills, baboons, and geladas) exhibit hallmarks of a high variance in male reproductive success, including a female-biased adult sex ratio and prominent sexual dimorphism. To explore the potential genomic consequences of such sex differences, we used a reduced representation genome sequencing approach to quantifying polymorphism at sites on autosomes and sex chromosomes of the tonkean macaque (Macaca tonkeana), a species endemic to the Indonesian island of Sulawesi. The ratio of nucleotide diversity of the X chromosome to that of the autosomes was less than the value (0.75) expected with a 1:1 sex ratio and no sex differences in the variance in reproductive success. However, the significance of this difference was dependent on which outgroup was used to standardize diversity levels. Using a new model that includes the effects of varying population size, sex differences in mutation rate between the autosomes and X chromosome, and GC-biased gene conversion (gBGC) or selection on GC content, we found that the maximum-likelihood estimate of the ratio of effective population size of the X chromosome to that of the autosomes was 0.68, which did not differ significantly from 0.75. We also found evidence for 1) a higher level of purifying selection on genic than nongenic regions, 2) gBGC or natural selection favoring increased GC content, 3) a dynamic demography characterized by population growth and contraction, 4) a higher mutation rate in males than females, and 5) a very low polymorphism level on the Y chromosome. These findings shed light on the population genomic consequences of sex differences in the variance in reproductive success, which appear to be modest in the tonkean macaque; they also suggest the occurrence of hitchhiking on the Y chromosome.

Lack of population genetic structure and host specificity in the bat fly, Cyclopodia horsfieldi, across species of Pteropus bats in Southeast Asia.

BACKGROUND: Population-level studies of parasites have the potential to elucidate patterns of host movement and cross-species interactions that are not evident from host genealogy alone. Bat flies are obligate and generally host-specific blood-feeding parasites of bats. Old-World flies in the family Nycteribiidae are entirely wingless and depend on their hosts for long-distance dispersal; their population genetics has been unstudied to date. METHODS: We collected a total of 125 bat flies from three Pteropus species (Pteropus vampyrus, P. hypomelanus, and P. lylei) from eight localities in Malaysia, Cambodia, and Vietnam. We identified specimens morphologically and then sequenced three mitochondrial DNA gene fragments (CoI, CoII, cytB; 1744 basepairs total) from a subset of 45 bat flies. We measured genetic diversity, molecular variance, and population genetic subdivision (FST), and used phylogenetic and haplotype network analyses to quantify parasite genetic structure across host species and localities. RESULTS: All flies were identified as Cyclopodia horsfieldi with the exception of two individuals of Eucampsipoda sundaica. Low levels of population genetic structure were detected between populations of Cyclopodia horsfieldi from across a wide geographic range (~1000 km), and tests for isolation by distance were rejected. AMOVA results support a lack of geographic and host-specific population structure, with molecular variance primarily partitioned within populations. Pairwise FST values from flies collected from island populations of Pteropus hypomelanus in East and West Peninsular Malaysia supported predictions based on previous studies of host genetic structure. CONCLUSIONS: The lack of population genetic structure and morphological variation observed in Cyclopodia horsfieldi is most likely due to frequent contact between flying fox species and subsequent high levels of parasite gene flow. Specifically, we suggest that Pteropus vampyrus may facilitate movement of bat flies between the three Pteropus species in the region. We demonstrate the utility of parasite genetics as an additional layer of information to measure host movement and interspecific host contact. These approaches may have wide implications for understanding zoonotic, epizootic, and enzootic disease dynamics. Bat flies may play a role as vectors of disease in bats, and their competence as vectors of bacterial and/or viral pathogens is in need of further investigation.

Scale-dependent effects of a heterogeneous landscape on genetic differentiation in the Central American squirrel monkey (Saimiri oerstedii).

Landscape genetic studies offer a fine-scale understanding of how habitat heterogeneity influences population genetic structure. We examined population genetic structure and conducted a landscape genetic analysis for the endangered Central American Squirrel Monkey (Saimiri oerstedii) that lives in the fragmented, human-modified habitats of the Central Pacific region of Costa Rica. We analyzed non-invasively collected fecal samples from 244 individuals from 14 groups for 16 microsatellite markers. We found two geographically separate genetic clusters in the Central Pacific region with evidence of recent gene flow among them. We also found significant differentiation among groups of S. o. citrinellus using pairwise F(ST) comparisons. These groups are in fragments of secondary forest separated by unsuitable "matrix" habitats such as cattle pasture, commercial African oil palm plantations, and human residential areas. We used an individual-based landscape genetic approach to measure spatial patterns of genetic variance while taking into account landscape heterogeneity. We found that large, commercial oil palm plantations represent moderate barriers to gene flow between populations, but cattle pastures, rivers, and residential areas do not. However, the influence of oil palm plantations on genetic variance was diminished when we restricted analyses to within population pairs, suggesting that their effect is scale-dependent and manifests during longer dispersal events among populations. We show that when landscape genetic methods are applied rigorously and at the right scale, they are sensitive enough to track population processes even in species with long, overlapping generations such as primates. Thus landscape genetic approaches are extremely valuable for the conservation management of a diverse array of endangered species in heterogeneous, human-modified habitats. Our results also stress the importance of explicitly considering the heterogeneity of matrix habitats in landscape genetic studies, instead of assuming that all matrix habitats have a uniform effect on population genetic processes.

Sex-linked inheritance in macaque monkeys: implications for effective population size and dispersal to Sulawesi.

Sex-specific differences in dispersal, survival, reproductive success, and natural selection differentially affect the effective population size (N(e)) of genomic regions with different modes of inheritance such as sex chromosomes and mitochondrial DNA. In papionin monkeys (macaques, baboons, geladas, mandrills, drills, and mangabeys), for example, these factors are expected to reduce N(e) of paternally inherited portions of the genome compared to maternally inherited portions. To explore this further, we quantified relative N(e) of autosomal DNA, X and Y chromosomes, and mitochondrial DNA using molecular polymorphism and divergence information from pigtail macaque monkeys (Macaca nemestrina). Consistent with demographic expectations, we found that N(e) of the Y is lower than expected from a Wright-Fisher idealized population with an equal proportion of males and females, whereas N(e) of mitochondrial DNA is higher. However, N(e) of 11 loci on the X chromosome was lower than expected, a finding that could be explained by pervasive hitchhiking effects on this chromosome. We evaluated the fit of these data to various models involving natural selection or sex-biased demography. Significant support was recovered for natural selection acting on the Y chromosome. A demographic model with a skewed sex ratio was more likely than one with sex-biased migration and explained the data about as well as an ideal model without sex-biased demography. We then incorporated these results into an evaluation of macaque divergence and migration on Borneo and Sulawesi islands. One X-linked locus was not monophyletic on Sulawesi, but multilocus data analyzed in a coalescent framework failed to reject a model without migration between these islands after both were colonized.

Range-wide mtDNA phylogeography yields insights into the origins of Asian elephants.

Recent phylogeographic studies of the endangered Asian elephant (Elephas maximus) reveal two highly divergent mitochondrial DNA (mtDNA) lineages, an elucidation of which is central to understanding the species's evolution. Previous explanations for the divergent clades include introgression of mtDNA haplotypes between ancestral species, allopatric divergence of the clades between Sri Lanka or the Sunda region and the mainland, historical trade of elephants, and retention of divergent lineages due to large population sizes. However, these studies lacked data from India and Myanmar, which host approximately 70 per cent of all extant Asian elephants. In this paper, we analyse mtDNA sequence data from 534 Asian elephants across the species's range to explain the current distribution of the two divergent clades. Based on phylogenetic reconstructions, estimates of times of origin of clades, probable ancestral areas of origin inferred from dispersal-vicariance analyses and the available fossil record, we believe both clades originated from Elephas hysudricus. This probably occurred allopatrically in different glacial refugia, the alpha clade in the Myanmar region and the beta clade possibly in southern India-Sri Lanka, 1.6-2.1Myr ago. Results from nested clade and dispersal-vicariance analyses indicate a subsequent isolation and independent diversification of the beta clade in both Sri Lanka and the Sunda region, followed by northward expansion of the clade. We also find more recent population expansions in both clades based on mismatch distributions. We therefore suggest a contraction-expansion scenario during severe climatic oscillations of the Quaternary, with range expansions from different refugia during warmer interglacials leading to the varying geographical overlaps of the two mtDNA clades. We also demonstrate that trade in Asian elephants has not substantially altered the species's mtDNA population genetic structure.

Resolution of the Hylobates phylogeny: congruence of mitochondrial D-loop sequences with molecular, behavioral, and morphological data sets.

Gibbons of the genus Hylobates likely speciated very rapidly following isolation by rising sea levels during the Pleistocene. We sequenced the hypervariable region I (HV-I) of the mitochondrial D-loop to reconstruct the phylogeny of this group. Although the results clearly supported monophyly of each of the six species, the relationships among them were not clearly resolved by these data alone. A homogeneity test against published data sets of a coding mitochondrial locus (ND3-ND4 region), behavioral characters (vocalizations), and morphological traits (including skeletal and soft tissue anatomy) revealed no significant incongruence, and combining them resulted in a phylogenetic tree with much stronger support. The Kloss's gibbon (H. klossii), long considered a primitive taxon based on morphology, shares many molecular and vocal characteristics with the Javan gibbon (H. moloch), and appear as the most recently derived species. The northernmost species (H. lar and H. pileatus) are the most basal taxa. These data suggest that ancestral gibbons radiated from north to south. Unlike other markers, the HV-I region can accurately identify members of different gibbon species much like a DNA barcode, with potential applications to conservation.

A complete species-level phylogeny of the Hylobatidae based on mitochondrial ND3-ND4 gene sequences.

The Hylobatidae (gibbons) are among the most endangered primates and their evolutionary history and systematics remain largely unresolved. We have investigated the species-level phylogenetic relationships among hylobatids using 1257 bases representing all species and an expanded data set of up to 2243 bases for select species from the mitochondrial ND3-ND4 region. Sequences were obtained from 34 individuals originating from all 12 recognized extant gibbon species. These data strongly support each of the four previously recognized clades or genera of gibbons, Nomascus, Bunopithecus, Symphalangus, and Hylobates, as monophyletic groups. Among these clades, there is some support for either Bunopithecus or Nomascus as the most basal, while in all analyses Hylobates appears to be the most recently derived. Within Nomascus, Nomascus sp. cf. nasutus is the most basal, followed by N. concolor, and then a clade of N. leucogenys and N. gabriellae. Within Hylobates, H. pileatus is the most basal, while H. moloch and H. klossii clearly, and H. agilis and H. muelleri likely form two more derived monophyletic clades. The segregation of H. klossii from other Hylobates species is not supported by this study. The present data are (1) consistent with the division of Hylobatidae into four distinct clades, (2) provide the first genetic evidence for all the species relationships within Nomascus, and (3) call for a revision of the current relationships among the species within Hylobates. We propose a phylogenetic tree as a working hypothesis against which intergeneric and interspecific relationships can be tested with additional genetic, morphological, and behavioral data.

Evolution of RAG-1 in polyploid clawed frogs.

Possible genetic fates of a gene duplicate are silencing, redundancy, subfunctionalization, or novel function. These different fates can be realized at the DNA, RNA, or protein level, and their genetic determinants are poorly understood. We explored molecular evolution of duplicated RAG-1 genes in African clawed frogs (Xenopus and Silurana) (1) to examine the fate of paralogs of this gene at the DNA level in terms of recombination, positive selection, and gene degeneration and in the absence of extensive recombination among alleles at different paralogs, (2) to test phylogenetic hypotheses about the origins of polyploid species. We found that recombination between different RAG-1 paralogs is infrequent, that degeneration of some paralogs has occurred via stop codons and frameshift mutations, and that this degeneration occurred in paralogs inherited from only one diploid progenitor species. Simulations and phylogenetic analyses of RAG-1 and mitochondrial DNA support one origin of extant tetraploids in Xenopus and at least one origin in Silurana, five allopolyploid origins of extant octoploids, and two allopolyploid origins of extant dodecaploids. In allopolyploid species, which inherit a complete genome from two different ancestors, genes inherited from the same ancestor have a longer period of coevolution than genes inherited from different ancestors. Because of this, gene ancestry could potentially influence gene fate: interacting paralogs derived from the same lower ploidy ancestor might have similar genetic destinies.

A mitochondrial DNA phylogeny of African clawed frogs: phylogeography and implications for polyploid evolution.

The African clawed frogs (Silurana and Xenopus), model organisms for scientific inquiry, are unusual in that allopolyploidization has occurred on multiple occasions, giving rise to tetraploid, octoploid, and dodecaploid species. To better understand their evolution, here we estimate a mitochondrial DNA phylogeny from all described and some undescribed species. We examine the timing and location of diversification, and test hypotheses concerning the frequency of polyploid speciation and taxonomy. Using a relaxed molecular clock, we estimate that extant clawed frog lineages originated well after the breakup of Gondwana, about 63.7 million years ago, with a 95% confidence interval from 50.4 to 81.3 million years ago. Silurana and two major lineages of Xenopus have overlapping distributions in sub-Saharan Africa, and dispersal-vicariance analysis suggests that clawed frogs originated in central and/or eastern equatorial Africa. Most or all extant species originated before the Pleistocene; recent rainforest refugia probably acted as "lifeboats" that preserved existing species, rather than "species pumps" where many new successful lineages originated. We estimate that polyploidization occurred at least six times in clawed frogs.